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It is well known that in a gas-filled duct or channel along which a temperature gradient is ap- 
plied, a thermal creep flow is created. Here we show that a mass and momentum flux can also be 
induced in a gas confined between two parallel structured surfaces at different temperatures, i.e. 
orthogonal to the temperature gradient. We use both analytical and numerical methods to compute 
the resulting fluxes. The momentum flux assumes its maximum value in the free-molecular flow 
regime, the (normalized) mass flux in the transition flow regime. The discovered phenomena could 
find applications in novel methods for energy-conversion and thermal pumping of gases. 

PACS numbers: 47.45.Dt free molecular flow, 44.15.+a,05.60.-k,51.10+y 



I. INTRODUCTION 

Thermally induced gas flow and Knudsen forces arising 
from the unbalanced momentum flux in a gas are among 
the most interesting phenomena in a rarefied gas with 
a non-uniform temperature field. The rapidly develop- 
ing micromachining technologies have made it possible 
to observe and apply these phenomena even at atmo- 
spheric conditions. For example, attractive and repulsive 
forces have been observed on heated microcantilevers in 
air PJ [5] , and a passive micromachined gas pump utiliz- 
ing thermal creep flow, entitled as the Knudsen pump, 
has been demonstrated [3J. The promising potential of 
these applications has called for a deeper understanding 
and study of transport phenomena in a rarefied gas. Sev- 
eral experimental and numerical studies have been con- 
ducted to study the Knudsen force on cantilevers @H2], 
radiometric forces [Ej , and Knudsen pumps [5] ■ 

Most Knudsen pumps proposed and developed up to 
now utilize thermal creep flow. The modeling of such 
flows dates back to 1879 when Maxwell jTUj derived a 
mathematical expression for the flow along a surface in- 
duced by a temperature gradient. Since then, the corre- 
sponding thermal creep flow through ducts or channels 
has been demonstrated in many experiments. The flow 
velocity is proportional to the temperature gradient, po- 
tentially requiring very large temperature spans when the 
channel length is no longer small, which has limited the 
application perspectives of thermal creep flow. Some at- 
tempts have been made to reduce the required tempera- 



ture span by considering a gas between parallel surfaces 
with a periodic, step-like topography [IX. Besides that, 
virtually no studies of transport phenomena in gases con- 
fined between two periodically structured surfaces seem 
to be available. 

It is well known that in a gas confined between two 
parallel isothermal walls of different temperatures no net 
mass and momentum fluxes parallel to the walls occur. 
However, as we will show in this article, with periodically 
structured surfaces it is possible to induce such fluxes. 
We consider a monoatomic gas governed by the Boltz- 
mann equation between two parallel surfaces equipped 
with periodic structures. The term "structure" has a 
twofold meaning. On the one hand, the surfaces can have 
a topography, on the other hand, there can be spatial 
modulations of the boundary condition expressing how 
a molecule gets reflected when impinging on a wall. We 
expect that especially when the Knudsen number Kn, 
being the ratio of the mean-free path of the molecules 
and the distance between the surfaces, is of the order 
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FIG. 1. (Color online) Schematic of the two-dimensional 
channel geometry considered made of two surfaces with tem- 
peratures Ti and T2. Different styles of the wall segments 
indicate different boundary conditions. On the right side, a 
corresponding geometry incorporating a wheel with suitable 
surface structures is sketched to suggest possible applications. 
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of one or larger, the surface topography and the spatial 
modulation of the boundary condition should have a sig- 
nificant influence on the phase space distribution of the 
molecules inside the channel and the resulting transport 
phenomena. We are mainly interested in the total mo- 
mentum and mass fluxes inside the gas that are created if 
the two surfaces are at different temperatures. The pur- 
pose of this work is to identify surface structures that are 
as simple as possible and nevertheless give rise to non- 
vanishing fluxes. The type of prototype geometry con- 
sidered is depicted in Fig. [T] together with the predicted 
directions of mass and momentum fluxes. The quanti- 
ties of interest are the total momentum and mass fluxes 
in a;-direction. It is clear that a non- vanishing flux along 
the x-direction will only result if the structure breaks the 
reflection symmetry with respect to x. As will become 
evident in the following, it suffices to structure one of the 
two surfaces. The top wall will therefore be considered 
flat and uniform. 



II. MOMENTUM FLUX 

Specifically, we are interested in the total momentum 
in x-direction per unit time the gas transfers between the 
top and the bottom surface. To further simplify matters, 
let us limit ourselves to the surface structures shown in 
Fig. [2] Besides choosing j3 — tt/2, the distribution of wall 
boundary conditions has been fixed. The red (solid) sec- 
tions denote a diffuse reflection boundary condition, the 
blue (dashed) ones a specular reflection boundary con- 
dition, whose mathematical form will be detailed below. 
Furthermore, for the time being we make the assump- 
tion that the mean-free path of the molecules is much 
larger than W, i.e. we assume gas dynamics in the free- 
molecular flow regime. We always suppose that the ther- 
mal conductivity of the wall material is so large that tem- 
perature variations along the surfaces can be neglected. 

Assuming no external force fields acting on the 
molecules, in the limit Kn — > 00 the solution of the 
Boltzmann equation can be reduced to tracing rays rep- 
resenting the paths of the molecules [T2Tll4j . We give a 
sketch of the derivation based on the approach presented 
in ref . [T5] . 

We denote by r% a point on the lower surface and 
by fi(r\,{-\) the distribution function of the incoming 
molecules incident with velocity ^ on T\. For those we 
have (£1 • Tii) < 0, where ri\ is inward normal at r\ (in 
this and the following the notation for the unit normals 
is analogous to the boundary points). The force exerted 
on the bottom surface by the incoming molecules is then: 



■ ni)/i(ri(si),£i)d 2 £idsi, 



£l-rn<0 



(1) 

where the integration is over an arc length parameter 
S\ through which the points at the bottom wall are de- 
fined via ri(si). Similarly, the recoil from the outgoing 
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FIG. 2. (Color online) The model geometry. Shown are the 
notations and an example of a tracing starting at 7*1 . 



molecules at 7*1 gives: 



' «i)/r(T , i(si),^i)d 2 ^idsi, 

£im>o 

(2) 

where /r(f"i>£i) is the phase space distribution of the 
outgoing molecules. The net force on the lower surface is 
then given by F — F t + F r . 

In the framework of the collisionless Boltzmann equa- 
tion /j(t"i,£i) is equal to / r (r,£i), where r is a point 
at the intersection of a straight line oriented along — £1 
with the boundary. The point r could be either on a 
diffusive or on a specular boundary, a situation depicted 
in Fig.[2]with the possibilities for r exemplified by rf, r d 
for the diffusive and by f for the specular wall segments, 
respectively. 

The boundary condition on the wall segments with dif- 
fuse reflection of molecules is expressed as 



f r (r d ,i 1 ) = v(r d )F 2D (r d ,i 1 ), 
2 / m 



F 2 > d ,a) 



/7T V '2T r d 



3/2 



exp 



mil 



(3) 



where m is the molecular mass, T r d the temperature of 
the wall at point r d (in energy units), and the num- 
ber of collisions from the incoming molecules at r d per 
unit length and unit time (the particle flux) v{r d ) has 
been factored out. This is a Maxwell boundary condi- 
tion [TU] with an accommodation coefficient of one. The 
particle flux is defined in the usual way as: v(r d ) = 

-f^. n «<o(z d - nd )Mr d ,i d )d 2 i d . 

From particle number conservation we have v(r d ) — 
/c 1 .n">o (& ■ nd ) fr{ rd ,ii)& 2 ii, from which the choice 
of the pre-exponential factor in Eq. Q is fixed. 

For the specular boundary segments the correspond- 
ing boundary condition is f r (f , £1) = fi(r, £1), where the 
expression on the right-hand side is obtained by further 
tracing in the direction of ~ii - The velocity £1 is related 
t° ii by reflection from the surface with a normal vector 
n, ii = ii — 2(n ■ ii)n. The specular boundary condi- 
tion is again Maxwell's boundary condition, but for a wall 
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with an accommodation coefficient of zero. Tracing par- 
ticle paths backwards from specular wall segments always 
brings us to points on a diffuse boundary, in our example 
in Fig. [2] the point r d at the upper wall, where the distri- 
bution function of reflected molecules is given by Eq. Q . 
By similar tracing procedures we find the relevant rela- 
tions for the distribution function in the equation for the 
recoil force Eq. ([2]). 

As a result, finding the phase space distribution of the 
molecules at points along the boundary of the domain 
essentially reduces to determining v(r). For illustration, 
a point r\ at the bottom wall receives particles traced 
to either a diffusive or to a specular wall segment. The 
particle flux's definition then gives an integral equation 

i/(ri) = - J (^■n 1 Hr d )F 2D (r d ,^)d 2 ^ 



n d (ri) 



(6-mMr d )F 2Z V,£i)d 2 £ 1 , (4) 

n(n) 

where i7 d (r 1 ) denotes the subset of incoming velocities 
£1 for which the tracing routine relates the point r\ to a 
point on a diffusive segment r d = r d (r i,£i), and simi- 
larly, 0(t"i) denotes the subset of incoming velocities for 
which the tracing routine relates the r\ to a point on the 
specular segment f, and eventually to a diffusive bound- 
ary point f d — f d (ri,^i). By direct substitution we find 
that this integral equation has a constant solution, i.e. 
v(r) = v, independent of position. 

To proceed, let us denote by Sq the projection of the tip 
of a fin onto the incline of its neighboring fin, cf. Fig. [2] 
For close-by fins or for sufficiently large values of a, so 
is positive, otherwise negative. This parameter allows to 
conveniently distinguish geometries such that a tracing 
procedure starting at ri and hitting the incline can only 
end at points on the top wall (sq < 0), and geometries for 
which the tracing after hitting the incline can end either 
on the top wall or on the vertical side wall of the lower 
surface (sq > 0). Evaluating the integrals in Eqs. 
and |2]) yields for one unit cell of length L : 



Fr 



where A 



Asm(a)y — — ajLcos(a), 
Asm(a) ^— — ctj L cos(a)- 



(*o > 0) 



-sq arctan 



•so 



L sin(a) 
j "V(a/5¥- VTT), and 



, (so < 0) 



[2m , r=- 



'To 



ttL 



F, 



tan(a) 



(5) 



(6) 



In the limit L 



00 we have F x oc 



sin(a) cos(a) H, as derived in [15]. As expected, the net 



force F x vanishes for equal temperatures T% = T2. Also, if 
the specular wall segments are replaced by diffuse bound- 
aries, F x vanishes. This shows that in the free-molecular 
flow regime a ^-momentum flux in the gas and the cor- 
responding force on the wall cannot be created by virtue 
of the wall topography alone. 

For comparison with the numerical calculations, we 
have considered the ratio F x /\F y \. Such a dimensionless 
ratio gives an immediate account of the expected magni- 
tude of the effect, since F y is the pressure force onto the 
wall. For example, when the free-molecular flow regime 
is approximated by considering a gas in a nanochanncl 
at standard conditions, the expression of Eq. ^ is the 
force resulting from a pressure of one atmosphere. 

To numerically solve the Boltzmann equation, Monte 
Carlo schemes were employed. Either a time-splitting 
Monte-Carlo method [El [17] with a hard-sphere col- 
lision model or a standard Direct Simulation Monte 
Carlo (DSMC) method with variable hard sphere col- 
lision model [18] was used. Both schemes utilize argon 
molecules as a gas species. 

DSMC is a particle method which tracks a number of 
simulated particles, each being a statistical representa- 
tion of a large cluster of real molecules. The molecular 
motions and intermolecular collisions are decoupled over 
a small time interval. Particles undergo a free convec- 
tion step followed by a collision step. The convection is 
treated deterministically while the intermolecular colli- 
sions are treated statistically. A detailed description of 
the DSMC method can be found in [15] . 

In all the simulations the cell size is set to be less than 
the mean free path and the time step is chosen so that 
no molecule can cross two cells within one time step. 
Molecules are initialized based on the average of the wall 
temperatures and the prescribed Knudsen number. A 
gas-wall interaction model to compute the post-collision 
molecular velocities as specified previously in this article 
is employed. 

In the case of the time-splitting method (the data in 
Figs. [3^a-b)), error bars were computed based on three 
independent runs with N ce u = 20, 40, 80 particles per 
computational cell. Following that, a Richardson ex- 
trapolation to N^ii — 00 was performed, allowing to 
estimate the error by computing the difference to the 
N ce u — 80 result. In the case of the standard DSMC 
method (the data in Figs. [3^c)-[5|, the error bars repre- 
sent the maximum deviations computed based on the last 
one third of total number of samples ranging from 68,000 
to 30,000,000. An average of N ce u — 40 particles per cell 
were used for these simulations. In all computations the 
temperatures were fixed to T t = 300 K, T 2 = 400 K. The 
consistency of the two numerical schemes was checked 
by comparing the results for F x for specific geometric 
parameter values at Kn = 1 and at Kn = 00. The data 
points obtained with the different methods were found 
to agree, with deviations within the error bars obtained 
by the extrapolation, and, at Kn = 00 , both were within 
half of their corresponding standard deviation from the 
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FIG. 3. (Color online) (a) and (b) Analytical (lines) and nu- 
merical (points) results for the force ratio as a function of geo- 
metric parameters. The circles represent the results obtained 
for Kn = co, the diamonds those for Kn = 1. (a) shows the 
dependence on a with L/W = 5/16, H/W = tan(a)/4. (b) 
displays the dependence on L/W, with H/W = 3/4, a = n/3. 
Inset: the data points at small values of L/W in more detail, 
(c) Numerical results for the force ratio as a function of Knud- 
scn number at a = tt/3, H/W = V3/4,L/W = 5/16. The red 
square indicates the prediction from the analytical model. 



analytical prediction, the red square in Fig. [3jc) . 

There are four geometric parameters, FL, L, W, a, and 
for the comparison between analytical and numerical cal- 
culations we have chosen two ways to fix them, as to have 
only one independent parameter to be varied. The first 
one is to fix the distance between the fins L and their 
base length H/ tan(a) and to change a, the second one 
is to vary L while keeping the other parameters the same. 
Figures [3] (a) and (b) show the analytical and numerical 
values of F x /F y obtained that way. Both Kn = 1 and 
Kn = oo were considered in the numerical studies. In 
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FIG. 4. (Color online) The area averaged s-velocity as a 
function of the Knudsen number. The averaging is done over 
a patch connecting the tip of a fin with the upper surface. 



addition, Fig. [3^c) shows the dependence of the force ra- 
tio on the Knudsen number. 

First of all it is worth noting that the analytical curves 
are well reproduced by the numerical data obtained for 
Kn = oo. When the Knudsen number is decreased to 
one, the data points still follow the same trend, but the 
magnitude of the dimensionless force is reduced by about 
30%. We therefore conclude that the analytical solution 
of the collisionless Boltzmann equation still makes quali- 
tatively correct predictions in the transition flow regime. 
The asymmetry of the fins (including the effect of the 
boundary condition) is the origin of the x-momentum 
flux inside the gas and the corresponding force onto the 
walls. Therefore it is expected that the maximum force 
values are obtained in a region around L/W — 1, and 
that the force decreases as the fins become sparser. This 
is exactly the behavior displayed in Fig. [3] (b). Fig. [3] (a) 
shows that the force is in fact significant: At Kn = 1 max- 
imum values of well above 1% are obtained. As expected, 
the force vanishes when the gas approaches a state where 
local thermal equilibrium prevails (cf. Fig. [31(c)). The 
maximum force values are reached in the free molecular 
flow regime. 

From an energetic point of view, when F x acts to dis- 
place the two parallel surfaces with respect to each other, 
heat is transformed into mechanical energy, as indicated 
in the schematic on the right-hand side of Fig. [I] We 
wish to emphasize that this is a novel energy conversion 
scheme that works without the usual volumetric expan- 
sion. Rather than functioning in a cycle, it allows to 
continuously extract mechanical work from a tempera- 
ture difference. 



III. MASS FLUX 

In the limit Kn — > oo where the force becomes maxi- 
mal, the mass flux in the gas vanishes. Such a scenario 
bears some analogy with results that have been obtained 
for closed domains [119 . It also shows that the force is 
not due to the shear stress caused by a fluid flowing over 
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FIG. 5. (Color online) The velocity field at Kn = 0.1. 
The background color (grayscale) contours indicate the cor- 
responding temperature field. 



a surface. At Kn = oo the phase space distribution 
of the gas is such that the expectation value of the x- 
momentum flux transferred between the walls is nonzero. 
The molecules emerging from the lower surface carry a 
negative x-momcntum flux, corresponding to the posi- 
tive value of F x the surface experiences. Away from the 
limiting case of infinite Knudsen number, some of these 
molecules collide with other gas molecules, thereby trans- 
ferring ^-momentum to the gas. It is expected that by 
such a mechanism the gas in the upper portion of the 
domain into which the fins do not extend is set into mo- 
tion. This picture suggests that the corresponding mass 
flux should be opposite to F x . 

One can also examine the picture from the other end 
of the gas regime. In the limit of Kn — > 0, the local 
equilibrium conditions ensure that both momentum and 
mass fluxes are zero and hence no force and mass flow 
exist. As Kn increases, non-equilibrium conditions start 
to prevail. Since the temperature gradient around the fin 
tips is likely to be the most significant, due to both the 
sharp tip and the asymmetric boundaries, one would ex- 
pect that deviations from local thermal equilibrium first 
occur in these regions, resulting in a larger non-zero mo- 
mentum flux compared to the gas in other regions. The 
flow induced by this momentum flux is then expected to 
have its maximum velocity near the tip and to gradually 
decrease towards the top wall, resulting in a flow pattern 
similar to Couette flow. As Kn increases further, the 
non-equilibrium region expands into the whole channel. 
The tip flow is no longer dominant and the flow pattern 



is dictated by the overall configuration of the channel. 

Figure [4] displays numerical results for the average flow 
velocity as a function of Kn. For these calculations a 
specific geometry with a = n/3,H/W = \/3/4, and 
L/W = 5/16 was chosen (the same geometry as for the 
next to last data point in Fig. |3ja)). As expected from 
the analytical results above and the Navier-Stokes limit 
of the Boltzmann equation, there is no flow in the limits 
Kn — > and Kn — > oo. The maximum average flow veloc- 
ity occurs in the region around Kn = 0.1. The behavior 
at large Knudsen numbers suggests that the mass flux 
can be interpreted as a perturbation effect to the zero 
mass flux, nonzero momentum flux phase space distri- 
bution derived in the previous section: The momentum 
transferred from particle-particle collisions sets the gas 
into motion. Figure [5] shows the flow pattern obtained at 
Kn = 0.1. The background color contours indicate the 
corresponding temperature field. The flow is in the nega- 
tive x-direction and a Couette-like flow pattern with the 
maximum velocity being near the fin tip is clearly visible, 
consistent with the simple picture drawn in the previous 
paragraph. Significant flow velocities of close to 1 m/s 
are reached, making the described principle relevant for 
the pumping of gases. 



IV. CONCLUSIONS AND OUTLOOK 

We have studied the transport processes in a gas be- 
tween a flat and a structured surface at different temper- 
atures. We have shown that with suitably chosen surface 
structures, both a momentum and a mass flux parallel to 
the surfaces can be induced. The momentum flux takes 
its maximum values in the free-molecular flow regime, 
the mass flux (in terms of the average velocity) in the 
transition flow regime. While it is a priori clear that 
suitable structures need to break the reflection symme- 
try, a purely geometrical symmetry breaking based on 
diffusively reflecting walls is not sufficient. The observed 
fluxes only occur because a suitable texture comprising 
surface patches with diffuse and specular reflection con- 
ditions have been selected. The described effects may 
find applications in micro-/nanoscale energy converters 
and micro-/nanoscale pumping technology. 
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